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When a fluid with a bulk liquid-vapor critical point is placed inside a static external field with spatial 
periodic oscillations in one direction, the bulk critical point splits into two new critical points and a triple 
point. This phenomenon is called laser-induced condensation [Mol. Phys. 101, 1651 (2003)], and it occurs 
when the wavelength of the field is sufficiently large. The critical points mark the end of two coexistence 
regions, namely between (1) a vapor and stackcd-fluid phase, and (2) a stacked-fluid and liquid phase. The 
stacked-fluid or "zebra" phase is characterized by large density oscillations along the field direction. We study 
the above phenomenon for a mixture of colloids and polymers using density functional theory and computer 
simulation. The theory predicts that the vapor-zebra and liquid-zebra surface tensions are extremely small. 
Most strikingly, however, is the theoretical finding that at their respective critical points, both tensions 
vanish, but not according to any critical power law. The solution to this apparent paradox is provided by the 
simulations. These show that the field divides the system into effectively two-dimensional slabs, stacked on 
top of each other along the field direction. Inside each slab, the system behaves as if it were two-dimensional, 
while in the field direction the system resembles a one-dimensional Ising chain. 

PACS numbers: 61.20.Gy,68.05.-n,82.70.Dd 
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I. INTRODUCTION 

Binary mixtures of sterically-stabilized colloids and 
non-adsorbing globular polymers are valuable model sys- 
tems. In particular, the addition of polymers causes 
an effective depletion attraction between the colloids. 
This can induce liquid-vapor type transitions in these 
systems, in much the same way as in an atomic 
fluid. Indeed, liquid-vapor demixing in colloid-polymer 
mixtures has been routinely studied in the past, us- 
ing theory, computer simulation, and experiments*— . 
More phenomena for which colloid-polymer mixtures 
are ideal model systems include equilibrium clustering 7 , 
"attractive" glasses^ 9 -, gelation^ ., numerous interfa- 
cial phenomena^— including capillary wave o 15 ' 16 , and 
wettingiZr— . 

Presumably the simplest model of a colloid-polymer 
mixture is the one proposed by Asakura and Oosawa 
(AO) 21-23 . Despite its simplicity, this model captures the 
essential physics 3 -! 5 -, yet remains simple enough to allow 
for both theoretical investigations (based, for instance, on 
liquid state and density functional theor y 24 i 25 ) , and com- 
puter simulation o 26 ! 27 . In agreement with experiments, 
the AO model features a bulk liquid- vapor critical point 2 , 
a freezing transition 2 ^, and also confinement effects at a 
single wall 29 ' 3 -, or between two parallel walls 3 -^— can be 
studied using this model. The AO model has also been 
used to study wetting2£~~, as well as phase separation 
in porous medial — . Also of interest is the phase be- 
havior of the AO model inside an external field, such as 
gravit y 46 ' 47 , or a spatially- varying periodic field 48 . 

In Ref. H|| a colloid-polymer mixture described by the 



AO model was subjected to a standing-wave external 
field, wavelength A, propagating in one direction (the in- 
direction in what follows). The field is thus effectively 
one-dimensional. This can be realized in experiments 
by placing a colloidal suspension inside a standing laser 
field 4 ^. While the influence of such a field on the freezing 
transition has been extensively studied, and is known 
to induce laser-induced freezin g 50 ' 51 ', relatively little is 
known regarding its effect on the liquid- vapor transition. 
Regarding the latter, Ref. HU proposes the following sce- 
nario: for sufficiently large A, the bulk liquid-vapor criti- 
cal point splits into two critical points and a triple point 
with an intermediate new phase which is partially con- 
densed in slabs perpendicular to the ^-direction. We call 
this phase the "zebra" phase in what follows. The results 
of Ref. HH were obtained by using fundamental measure 
density functional theory for a colloid-polymer mixture 
described by the AO model. 

In this work, we use the same model and technique to 
calculate the surface tensions between all three coexist- 
ing phases. We find that the vapor-zebra and liquid-zebra 
surface tensions are extremely small. Moreover, to our 
surprise, upon approach of the critical points, the latter 
tensions do not yield the expected critical power law be- 
havior. To clarify the nature of the critical points, we use 
Monte Carlo simulations and finite-size scaling. The sim- 
ulations confirm all the trends predicted by the theory, 
and also explain why the vapor-zebra and liquid-zebra 
surface tensions do not become critical. The main find- 
ing is that the standing-wave external field divides the 
system into effectively two-dimensional slabs, stacked on 
top of each other along the z-direction. Inside each slab, 
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the system behaves as if it were two-dimensional, while in 
the ^-direction the system resembles an effectively one- 
dimensional system. Hence, by turning on the external 
field, the bulk system is "split-up" into 2 + 1 separate 
dimensions. With this picture in mind, most of our ob- 
servations can be intuitively understood. 



II. MODEL AND UNIT CONVENTIONS 

To describe the interactions between colloids (c) 
and polymers (p) we use the Asakura-Oosawa (AO) 
model2i~— . The colloids and the polymer coils are both 
assumed to be spherical objects with respective diameters 
a and a p . In what follows, the colloid diameter a = 1 will 
be the unit of length, and the colloid-to-polymer size ra- 
tio is denoted q = a p / 'a. The interaction between colloid- 
colloid and colloid-polymer pairs is hard-core, while the 
polymer-polymer interaction is ideal, leading to the fol- 
lowing pair potentials 



U cc (r) 



*'C p 



(r) 



oo r < a 

otherwise, 

oo r < (a + <7 p )/2 

otherwise, 



(1) 



(2) 



u pp (r) = 0, (3) 

with r the center-to-center distance. As the interactions 
are either hard-core or ideal, the temperature T can be 
scaled out (it only sets the energy scale fc^T, where fc^ is 
the Boltzmann constant). We mostly use a grand canoni- 
cal ensemble, i.e. the system volume V, the colloid chem- 
ical potential p, c , and the polymer chemical potential fi p 
are fixed, but the number of colloids N c and polymers N p 
inside V fluctuates. The particle densities are defined as 
Pi = Ni/V, with i e (c,p). We also introduce the col- 
loid packing fraction rj c — p c v c , with v c = ira 3 /6 the 
volume of a single colloid. Following convention 2 -, we do 
not use the polymer chemical potential itself, but rather 
the polymer reservoir packing fraction rf p . It is defined 
as the packing fraction rj p = p p v p , v p = na p /Q : of a pure 
polymer system at given fi p (for the AO model, such a 
system is simply an ideal gas, and hence rf p oc e tJ/p ^ kBT ). 
In order to compare p c between theory and simulation, 
the thermal wavelength A = a j 2 in what follows. 

In the AO model, there is an effective attraction be- 
tween the colloids. This can be shown formally by "inte- 
grating out" the polymer o 21 ! 52 , which leads to an effec- 
tive colloid-colloid pair potential with an attractive well; 
the well-depth is proportional to rf p . Hence, the bulk AO 
model undergoes a liquid- vapor type transition, with rf p 
playing the role of inverse temperature (by bulk we ex- 
plicitly mean a three-dimensional system in the absence 
of any surfaces or fields) . 

The extension of this work is to consider the AO model 
inside an external field propagating along the z-direction 



with Vq the field amplitude, and A the wavelength. 
We emphasize that the field is static: it oscillates in 
space, not in time. Physically, Eq. ^ resembles a one- 
dimensional standing optical wave, which could easily be 
realized experimentally using a laser beam. In this work, 
we assume that Eq. ([J| acts on the colloidal particles 
only, but that the polymers do not "feel" the field. The 
Hamiltonian of the system is thus defined by the AO 
pair potentials, Eqs.(l)-(3), plus an external field contri- 
bution J2i=i Vcxt(zi), where the sum is over all colloids, 
and with z\ the z-coordinate of the i-th colloidal particle. 



III. DENSITY FUNCTIONAL THEORY 

In density functional theory (DFT), the equilibrium 
density profiles are the ones that minimize the grand 
canonical free energy functional 



Q[T, p c ,rj p ,p c (z),p p (z)}, 



(5) 



Vextiz) = V COS (2ttz/A) . 



(4) 



where p c {z) is the average colloid density at position z 
along the field direction, and p p (z) that of the polymers 
(we thus use an effectively one-dimensional set-up in our 
DFT calculations). In what follows, we also use the col- 
loid packing fraction profile T) c (z), which is simply p c (z) 
multiplied by v c . The system is symmetric around z = 0, 
and we include periodic boundary conditions. Based on 
the proof that a free energy functional SI indeed exists^, 
we use the fundamental measure approach of Ref. Hi] to 
approximate Eq. (JSJ) ; see Appendix |A"1 for full details. We 
numerically solve the resulting stationarity equation us- 
ing a Picard iteration scheme^. 



A. Phase diagram 

First, we revisit Ref. l48l and hence choose the size ratio 
q = 0.6, the wavelength of the external field A — 8.192, 
and its amplitude Vq = 0.5fcsT. Fig. Q] shows the phase 
diagram obtained from our DFT calculation. Here, we 
use the grand canonical representation, i.e. we plot the 
binodals in the (fx c ,r] p ) plane. Clearly visible is the 
characteristic "inverted letter Y" or "pitchfork" topology 
(solid lines). For comparison, the dashed line shows the 
bulk binodal, i.e. obtained without the external potential 
of Eq. ((4]). We first note that the bulk critical point oc- 
curs at a value of rf v below that of the vapor-zebra and 
liquid-zebra critical points. This is to be expected as con- 
finement generally lowers transition temperatures. The 
more striking feature of Fig. [T] is that r] p of the liquid- 
zebra critical point exceeds that of the vapor-zebra one: 
we find ?7p'™ » 0.526 and 77^. w 0.540. In contrast, in 
Ref. |H, no difference could be detected, which demon- 
strates the improved accuracy of the present work. The 
fact that the critical "inverse temperatures" differ is a 
genuine feature, since our simulations reveal the same 
trend. As rj p increases, the binodals approach each other, 
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FIG. 1. Phase diagram of the AO model with q — 0.6 in 
(/u c ,7?p) representation as obtained using DFT. The dashed 
curve shows the bulk binodal, the solid curves the binodals 
inside the external field of Eq. @ using A = 8.192, and 
Vo = 0.5fcsT. The lower three dots indicate critical points, 
the upper dot the triple point. The open squares denote the 
state points at which the colloid density profiles of Fig. were 
measured. 



i 1 1 1 r 




FIG. 2. Same as Fig. [TJ but using q = 1.0, with external field 
parameters A = 10, and Vo = 0.4 ksT . 



and at the triple point, rf v tr « 0.755, they meet. In order 
to facilitate the comparison to computer simulation later 
on, we also present the phase diagram for q = 1, using 
field parameters A = 10, and V = 0Ak B T (Fig. [3]). We 
obtain the same overall topology, but the region where 
the zebra phase occurs has broadened. In addition, the 
binodals are shifted to significantly lower colloid chemical 
potential. 

Next, we consider the structural properties of the 
phases. The key difference with the bulk AO model is 
that, in addition to a vapor and liquid phase, we now also 
have the zebra phase. The latter phase arises when rj^ 
is chosen between the critical and triple points, and with 
the colloid chemical potential chosen appropriately. To 
characterize the phases, we have measured colloid den- 
sity profiles rj c (z) along the direction of the laser field at 
three points in the phase diagram of Fig. [TJ indicated by 




FIG. 3. Equilibrium colloid density profiles r] c (z) measured 
along the direction of the laser field for the three state points 
indicated by open squares in the phase diagram of Fig. [TJ 
corresponding to the vapor phase (fic/ksT — 6.72), the zebra 
phase ((Jc/fceT = 7.12), and the liquid phase (n c /kBT = 
7.52). The profiles were obtained at fixed r)p = 0.65. 

open squares. The latter correspond, from left to right, 
to a vapor state, a zebra state, and a liquid state. The 
density profiles are shown in Fig. [3J The salient feature 
is that all three phases display density modulations in 
the ^-direction, but the average density and amplitude 
differ. The average colloid density is low in the vapor 
phase, high in the liquid phase, with only a modest den- 
sity amplitude in both phases. The most striking feature 
of the zebra phase is the unusually large density ampli- 
tude, which oscillates between the average density of the 
vapor and liquid phase. 

B. Interfaces and interfacial free energies 

We now consider the interfaces between the coexisting 
phases and the corresponding surface tensions. Above 
the triple point, rj^ > 77^ tr , liquid and vapor coexist, with 
a corresponding liquid- vapor surface tension 7i v . To cal- 
culate 7i v , we first compute the equilibrium colloid den- 
sity profiles r) c {z) of the pure vapor and liquid phase; the 
latter yield the Gibbs free energies ^v,pure and ^i,pure: 
respectively (at coexistence: f2 v ,pure = £\purc = ^pure). 
The density profiles of the pure phases schematically re- 
semble those of Fig. [3J Next, we consider J] c (z) of a sys- 
tem containing a liquid-vapor interface, from which we 
obtain Jlivjnt- An example is shown in Fig. |D(a), using 
two values of 77^. For sufficiently large ?7p, we observe 
small oscillations at high densities close to the interface. 
Since the surface tension is the excess free energy per 
unit of area, it follows that 

^purc ^lv,int / n \ 

7iv = u , (6) 

where A is the area of the interface (the factor 1/2 results 
from the fact that two interfaces are present in our DFT 
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FIG. 5. Variation of the surface tensions 7i v , 7 VZ , and 7i z 
with rf p . Note in particular the extremely small values of 7 VZ 
and 7i z . These data were obtained using our DFT for the 
AO model with q = 1, and external field parameters A = 10, 
V = 0.4 k B T. 



to those obtained without the interface! From this obser- 
vation, and comparing to Eq. ([5]) , one can already deduce 
that 7vz,7iz must be very small. 

In Fig. [51 we summarize the results of the DFT surface 
tension calculations, where the various tensions are plot- 
ted as function of rf^. Starting above the triple point, yi v 
is finite; by decreasing rf p , j\ v vanishes at the triple point. 
At the triple point, 7 VZ and 7i z are finite; by decreasing if p 
further, 7 VZ vanishes at rf p — f/p ^ of the vapor-zebra crit- 
ical point, while 7i z vanishes at 77^ = j^'k of the liquid- 
zebra critical point. Note the extremely small values of 
7 VZ and 7i z over the entire range between the critical and 
triple points. In fact, it always holds that 7i v > 7 VZ + 7iz, 
which implies there is no complete wetting of the "zebra" 
phase for rf p > rf p >tr . 



Critical behavior 



FIG. 4. Equilibrium colloid density profiles r) c (z) showing the 
various interfaces for the AO model with q = 1, inside the 
external field of Eq. using field parameters A = 10, and 
Vo = 0.4fcsT. Shown is the liquid-vapor interface (a), the 
vapor-zebra interface (b), and the liquid-zebra interface (c), 
each time for two values of rfp as indicated. 



set-up). The calculation of the vapor-zebra surface ten- 
sion 7 VZ , and of the liquid-zebra surface tension 7i z , which 
become defined below the triple point, is performed anal- 
ogously. To this end, one needs to compute r/ c (z) for a 
system containing a vapor-zebra and liquid-zebra inter- 
face. Some typical density profiles are shown in Fig.[4jb) 
and (c). It is striking that the interfaces are extremely 
sharp: even very close to the interface position, the den- 
sity profiles of the coexisting phases are almost identical 



We now discuss the critical behavior of our equilibrium 
density profiles close to the vapor-zebra and liquid-zebra 
critical points. Since our DFT is a mean-field theory, we 
should recover mean-field critical exponents. This is not 
to suggest that the universality class of the AO model 
is the mean-field one - it is noti2& - but rather that we 
wish to test the internal consistency of our theory. To 
this end, we introduce the vapor-zebra order parameter 

A vz = — I r) CiZ (z) -rj CtV (z)dz, (7) 
A Jo 

where ?y CiZ (f? c ,v) denotes the equilibrium colloid density 
profile of the zebra (vapor) phase. Note that A vz above 
is just the difference between the average colloid density 
of the vapor and zebra phase. The liquid-zebra order 
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FIG. 6. Order parameters A vz and Ai z as function of the 
distance t from their respective critical points, where rfc™ ~ 
0.6817 and Vp'% ~ 0.717 were used. Note the double logarith- 
mic scales! The straight line corresponds to a power law with 
critical exponent ft = 1/2 of mean-field theory. For clarity, 
the data for Ai z have been shifted upward by half a decade. 
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FIG. 7. Surface tensions 7 VZ , 7i z , and the bulk liquid- vapor 
surface tension, as function of the distance t from the respec- 
tive critical points, where for Vp™ and r]p^ r the values of Fig. [6] 
were used, and fjp'cr ~ 0.6385. Note again the double loga- 
rithmic scales! The straight line corresponds to a power law 
with critical exponent /i = 3/2 of mean- field theory. The key 
message is that 7 VZ and 7i z do not conform to the mean-field 
critical exponent (we instead find an exponent fi w 0.9). 

parameter is defined analogously 

Aiz = ^ J Vc,i(z) - r) CtZ (z) dz. (8) 

Near the critical points, we expect power law decay of 
the order parameter 

A x oc^, i = r£ - r?^ > 0, xe(vz,lz), (9) 

with critical exponent (3 — 1/2 for mean- field theory. We 
compute these order parameters as function of r] p and 
plot them on double logarithmic scales in Fig.[6j where on 
the horizontal axes the distance from the critical point t 



is shown. The power law of Eq. ([9]) with mean-field ex- 
ponent /3 = 1/2 is strikingly confirmed! 

Next, we consider the critical behavior of the vapor- 
zebra and liquid-zebra surface tension 

7xOC^, xe(vz,lz), (10) 

with t > defined as above, and where the critical ex- 
ponent /i = 3/2 for mean- field theory. In Fig. [7J we 
plot both surface tensions as function of t, again using 
double logarithmic scales (this plot is simply a rescaling 
of the data of Fig. [5|. For completeness, we also show 
the liquid-vapor surface tension of the bulk AO model, 
i.e. in the absence of the laser field. The puzzling result 
is that, while the bulk tension conforms to /i = 3/2 as 
expected, we do not recover the expected mean-field criti- 
cal exponent for the vapor-zebra and liquid-zebra surface 
tensions. From this we conclude that 7 VZ and 7i z do not 
become critical. We postulate there must be a "hidden" 
surface tension 7h which instead conforms to Eq. (|T0|) at 
the vapor-zebra and liquid-zebra critical points; finding 
the corresponding "hidden" surface is one of the chal- 
lenges facing the simulations. 

IV. MONTE CARLO RESULTS 

We now use computer simulations to corroborate the 
DFT findings, and to shed light on the peculiar nature 
of the vapor-zebra and liquid-zebra critical points. We 
simulate the AO model (defined in Section [TTJ) inside the 
external potential of Eq. (0} using grand canonical Monte 
Carlo^. In the grand canonical ensemble, the colloid 
chemical potential /i c and the polymer "reservoir packing 
fraction" rf v are fixed, while the number of colloids N c 
and polymers N p in the system fluctuate. We remind 
the reader that r£ plays the role of inverse temperature. 
To simulate efficiently, a grand canonical cluster move is 
useo~, combined with a biased sampling scheme^. The 
simulations are performed in a V — L x x L y x L z box with 
periodic boundaries. The laser field, Eq. ((3]), propagates 
along the edge L z of the box, and hence we choose L z — 
nX, with integer n > 0, and A the wavelength of the field. 
In what follows, the colloid-to-polymer size ratio q = 1, 
A = 10, and the laser field amplitude Vo/fc^T = 0.4. 
The key output of the simulations is the order parameter 
distribution P(r] c ) (OPD) defined as the probability to 
observe the system in a state with colloid packing fraction 
i] c . From the (normalized) OPD, one readily computes 
the average colloid packing fraction 

(Vc) = J VoP(Vo)dT)o, (11) 

as well as the colloidal compressibility 

Xc = V-((r, c 2 )-(^ c ) 2 ), (12) 
and the Binder cumulant££ 

Q= (m 2 ) 2 /(™ 4 ), m = T? c -<Tfc>. (13) 
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FIG. 8. This figure illustrates how the vapor-zebra (VZ) and 
liquid-zebra (LZ) transitions are located in grand canonical 
simulations; rf v = 1.0, L x = L y = 12, and L z = 2A = 20 
are used in both plots, (a) The colloidal compressibility Xc 
versus the colloid chemical potential /i c ; the chemical poten- 
tial of the left (right) peak yields ). (b) The cumu- 
lant Q as function of (rjc) (solid curve). The dashed curve 
shows the compressibility \c on an arbitrary vertical scale. 
The compressibility maxima coincide with maxima in Q, the 
adjacent minima of which are labeled A, B ,C, D. The hori- 
zontal double-arrows mark the vapor-zebra and liquid-zebra 
order parameters. 
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FIG. 9. Phase diagram of the AO model inside the laser field 
of Eq. (Q in grand canonical representation (a) and reservoir 
representation (b); dots mark the critical points obtained via 
finite-size scaling. No scaling analysis was performed to locate 
the triple point, but based on the merging of the curves in (a), 
and also on the shape of the OPD (Fig. I17[) . we expect that 
^p.tr ~ 1-21 — 1-22 (horizontal lines). The dashed curves in 
(b) are power laws corresponding to the critical exponent ft 
obtained in Fig. UOf b). As the system size is increased, the 
binodals obtained in the finite system (solid curves) smoothly 
approach these power laws. Also labeled in (b) are the various 
coexistence regions. 



We emphasize that the above quantities, as well as the 
OPD, depend on all the model parameters, in particular 
the system size, the imposed colloid chemical potential 
fi c , and the "inverse temperature" rf v . 

A. Phase diagram 

To obtain the phase diagram, we vary the colloid chem- 
ical potential /i c at fixed "inverse temperature" rf p ; phase 
transitions correspond to peaks in the colloidal compress- 
ibility. An example is shown in Fig.[5ta), where \c versus 
fi c is plotted. We observe two sharp peaks, indicating 
two transitions. The left (right) peak corresponds to the 
vapor-zebra (liquid-zebra) transition, and from the peak 
position /i™ (/ijf) can be "read-off". Note that, above 



the triple point (rjl > ?7p tr ), Xc versus /i c reveals only 
one peak, then corresponding to a liquid-vapor transi- 
tion. For a range of rf vl we record the value (s) of the 
colloid chemical potential where \c is maximal, and plot 
these as points in the (/z c , 77^) plane. The resulting phase 
diagram is shown in Fig.^a), and the "inverted letter Y" 
topology predicted by the DFT is strikingly confirmed. 

The dots in Fig. [9] mark the vapor-zebra and liquid- 
zebra critical points, which we obtained using finite size 
scaling. For a given value of Tyjj and system size, we vary 
the colloid chemical potential fi c , and record the average 
colloid packing fraction (r] c ), the colloidal compressibility 
Xc, and the Binder cumulant Q. A typical result is shown 
in Fig.[5Jb), where Q and Xc versus (ry c ) are plotted (these 
curves are thus parametrized by /i c ). The key message of 
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FIG. 10. Finite-size scaling analysis of the vapor-zebra critical 
point. In both plots, L x — L y = L is varied, while L z — 2A 
is fixed, (a) The cumulant Q vz versus rf p for different L; the 
intersection yields J/p™. Note that, to the left (right) of the 
intersection, the cumulants approach 1/3 (0), in agreement 
with Eq. ()16[) . (b) Scaling plot of the vapor-zebra order pa- 
rameter. Plotted is AvzL' 3 /" versus tL 1 ^ , where v and /3 were 
tuned until a good collapse of the data was observed. 



Fig. Hlb) is that the compressibility maxima of the vapor- 
zebra and liquid-zebra transitions coincide with maxima 
in the cumulant. Adjacent to the cumulant maximum 
of the vapor-zebra transition, we observe two minima, 
indicated by the points 

A = «77c)7z, Q7z) > B ee «r /c )+ , Q+ ) , (14) 

and, similarly, for the liquid-zebra transition 

Cee(( Vc )- z ,Q- z ), D=((rk)+,Q+). (15) 

In the thermodynamic limit, it holds tha t 58 ' 59 



lim Q v 

L x jLy , L z — )> oo 



1/3 rf v <rf^ 
rf v > 



(16) 



with Q vz = (Q- + Q+ )/2, and the value of Vp 
at the vapor-zebra critical point. Hence, by plotting Q vz 
versus r) p for a number of different system sizes, curves for 



different system sizes intersect at r) p = T) p '™, which can 
be used to locate the critical point. Of course, to locate 
the liquid-zebra critical point, one analogously analyzes 

Qi» = (q; + q£)/2._ 

To perform the finite-size scaling analysis, we vary the 
lateral box extensions L x = L y = L, but keep the elon- 
gated extension fixed at L z = 2 A = 20. We assume that 
the divergence of the correlation length is "cut-off" in 
the z-direction by the laser field, and so we do not need 
to scale in this direction (this assumption will be justi- 
fied in the next section where the static structure factor 
is discussed). Since the correlations diverge only in the 
two lateral directions, the critical behavior is effectively 
two-dimensional. In Fig. ITUT a). we plot Q vz versus r) p 
for three values of L. In agreement with Eq. (I16|) . an 
intersection point is observed, from which we conclude 
that ?7p'™ ~ 0.975. A similar analysis of the liquid-zebra 
transition yields 77!'™. « 1.043 (not shown). It is striking 
that the scaling analysis confirms the DFT prediction 
Vp% > ^p.cr- To estimate the colloid chemical poten- 
tial /ic Z cr °f the vapor-zebra critical point in the thermo- 
dynamic limit, we measured /xjf 2 of the compressibility 
maximum at rf p — Tip'™ for finite L, and extrapolated to 
L — » 00 assuming /i™ cr — /i™ oc 1/L. In this extrapola- 
tion, we ignore all subtleties concerning field and pressure 
mixing^, but emphasize that such effects are tiny on the 
scale of the phase diagram in Fig. HJa). The resulting 
estimate reads as (j% z CI « 1.29, while for the liquid-zebra 
transition /j[. z cr m 1.75 is obtained. 

Next, we consider the scaling of the order parameter. 
The cumulant minima A and B of Fig. [5]Jb) readily yield 
A vz = (?7c}vz — (v^vz as order parameter for the vapor- 
zebra transition. In the vicinity of the critical point 
A vz oc r/ 3 , with t = (r? p - t££)/>#,Si t > 0, and crit- 
ical exponent j3. The result is shown in Fig. llOf b). where 
we used the standard finite-size scaling practice 61 of plot- 
ting A vz L /3 / I/ versus tL 1 /" , with v the correlation length 
critical exponent. Provided suitable values of 77^™, (3 and 
v are used, the data for different L collapse. Reasonable 
collapses can indeed be realized, using for 77^™ the cumu- 
lant intersection estimate of Fig.fTDTa). fijv ~ 0.25 — 0.35, 
and v ~ 0.85 — 1.10. An analysis of Ai z , which one ob- 
tains from the minima C and D of Fig.fSJb), yields similar 
results (not shown). 

Finally, with the critical point parameters known, it 
becomes possible to calculate the phase diagram in reser- 
voir representation, as is commonly done for the AO 
model. To this end, we record the colloid packing fraction 
of each of the cumulant minima A, B, C, D in Fig. [S^b) 
as function of i] p ; the latter "trace-out" a curve (bin- 
odal) in the (?7c,?7p) plane. The result is shown in 
Fig. (Hb), where dots again mark the critical points. 
To estimate the colloid packing fraction of the vapor- 
zebra critical point, we measured the finite-size "diam- 
eter" 5 L , VZ ee ((t7 c )7 z + (?7 c )+)/2, using if p = rfy™ and 
the colloid chemical potential /i™ of the compressibility 
maximum; the diameter was then extrapolated to L — > 00 
assuming <5oo iVZ — 5l,vz oc 1/L, which again ignores field 
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FIG. 11. The colloid-colloid static structure factors obtained 
at the vapor-zebra critical point, (a) The structure factor 
S\\ (q) measured in the direction of the laser field for several 
values of L z , and fixed L x = L y — 10. There is no divergence 
as q — > 0, indicating that critical fluctuations in the direction 
of the field are "cut-off", (b) The structure factor S±, z (q) 
measured in directions perpendicular to the laser field for sev- 
eral values of z; the system size used equals L x — L y = 20, 
and L z — 4A. The key message is that S± iZ (q) diverges as 
q — s> 0, but only for selected values of z. 



and pressure mixing effects^. In this way Soo.vz » 0.055 
is obtained, while an analogous procedure for the liquid- 
zebra critical point yields <5oo,i z « 0.203. 



B. Nature of the critical point 

A key assumption in the finite-size scaling analysis of 
the previous section is that the critical correlations are 
"cut-off" in the z-direction, i.e. the direction along which 
the laser field of Eq. Q propagates. To justify this as- 
sumption, we consider the colloid-colloid static structure 
factor S(q) — (^-| Y^=i exp(zg* ■ rj-) | 2 )i w i tn (') a thermal 
average, the sum over all j = 1, . . . , n c colloidal particles 
whose centers are inside a test volume v, and r*j the posi- 
tion of the j-th colloid. As usual, wavevectors are given 
by q = 2ir(k/L x , l/L y , m/L z ), integers k,l,m > 0, with 
the constraint that k + 1 + m ^ 0. We also introduce the 



wavevector magnitude q 2 — q- q. 

To probe the correlations in the z-direction, we use as 
test volume v a narrow cylinder, with a diameter equal to 
the colloid diameter, placed parallel to the z-axis; due to 
the symmetry of the system, the location where the cylin- 
der intersects the iry-plane is irrelevant. We then calcu- 
late the structure factor S\\(q), which is obtained using 
only the wavevectors q« = 2ir(0, 0, m/L z ). In Fig. ITTT a). 
we plot S\\ (q) measured at the vapor-zebra critical point. 
These data were obtained in a semi-grand canonical en- 
semble: the colloid packing fraction and ryjj are fixed to 
their critical values (?/ c = Soo. vz , = rf^ v £), while the 
number of polymers fluctuates. The important message 
to take from Fig. [TTTaj is that, in the limit g^fl, there 
is no sign of a divergence Sy^) — > oo. Hence, there are 
no critical fluctuations in the z-direction. Note that the 
peak at q « 0.63 corresponds precisely to 2n/X of the 
laser field. The analysis of Su (q) at the liquid-zebra crit- 
ical point leads to similar conclusions (not shown). 

Next, we consider the static structure factor S± z (q) 
measured in the lateral ^-directions, i.e. perpendicular 
to the laser field. In this case, wavevectors take the form 
q± = 2ir{k/L x , l/L v , 0), and as test volume v we use a 
narrow L x x L y x Az slab, placed parallel to the xy-plane 
at "height" z (the slab thickness Az equals one colloid di- 
ameter) . Since the system is not translation invariant in 
the z-direction, it matters at which z-coordinate the slab 
is located, and so S±, z (q) depends on z. In Fig.lTTTb). we 
plot S±^ z (q) at the vapor-zebra critical point for several z, 
again obtained using the semi-grand canonical ensemble. 
The key message to take from Fig. ITTTb) is that S± tZ (q) 
does diverge as q 0, but only for certain values of z. 
The analysis of S±. z (q) at the liquid-zebra critical point 
leads to similar conclusions (not shown). 

To summarize: From the static structure factor S\\(q), 
we conclude that critical fluctuations in the z-direction 
are absent. This justifies our previous assumption that 
the critical behavior is effectively two-dimensional, such 
that finite-size scaling may be performed by varying only 
the lateral box extensions L x = L y = L, while keep- 
ing L z fixed. The analysis of S± tZ (q) reveals that criti- 
cal fluctuations indeed develop in the lateral directions, 
but only at certain z values. The critical behavior is 
thus localized in effectively two-dimensional slabs per- 
pendicular to the laser field, "sandwiched" between slabs 
where the system is non-critical. To make this explicit, 
we show in Fig. I12f a) the variation of S± lZ (q m i n ) with 
z for the vapor-zebra critical point, where g m ; n = 2ir/L 
denotes the magnitude of the smallest accessible lateral 
wavevector. The figure strikingly shows that S± tZ (q m m) 
diverges with L only at selected z values, i.e. the criti- 
cal behavior is indeed spatially localized in slabs. Note 
that the critical slabs correspond to regions of enhanced 
colloid density: S±^ z (q m [ n ) is "in-phase" with the col- 
loid density profile i] c (z) (Fig. [12T b)). Interestingly, 
at the liquid-zebra critical point, this trend is reversed 
(Fig. USD- s ) cx Xc — , with Xc the col- 

loidal compressibility, we expect S± tZ (q m i n ) oc L 7 /" in 
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FIG. 12. Profiled quantities obtained at the vapor-zebra crit- 
ical point, (a) The variation of S± tZ (q m i M ) with z for several 
L. (b) The average colloid packing fraction rj c (z) measured 
along the z-direction (solid curve). The dashed curve shows 
the external laser potential of Eq. (j4| on an arbitrary vertical 
scale; regions dense in colloids coincide with minima of the 
potential. 



the critical slabs^i. Here, 7 is the compressibility criti- 
cal exponent; by fitting the peak values in Fig. Il"27 a) to 
this scaling law, 7/1/ ~ 1.3 — 1.4 is obtained. For the 
liquid-zebra critical point, a similar ratio is found, see 
Fig. [TBTa). where S± t z(q m in)L~ 7 ' L ' versus z is shown (in 
this scaled representation, the peak values for different L 
collapse) . It is reassuring that the critical exponent ratios 
obtained in our scaling analysis conform to hyperscaling, 
2/3/^ + 7/z/ = d = 2, as the reader can verify. Inter- 
estingly, our critical exponent ratios are rather different 
from 2D Ising values {fijv = 1/8, 7/1/ = 7/4), which we 
would naively have come to expect (only our v estimate is 
somewhat consistent with v = 1 of the 2D Ising model) . 



C. The coexistence region 

We now consider the vapor-zebra and liquid-zebra two- 
phase coexistence regions, see Fig. [9jb) , where the cor- 
responding transitions are first-order. To this end, we 
choose ?7p above the critical points, but still below the 



FIG. 13. The analogue of Fig. [12] for the liquid-zebra critical 
point. In (a), 7/1/ = 1.4 is used. 



triple point, and measure the OPD P(rj c ). In Fig. [TH we 
show P(t]c) using /i c = fj% z of the vapor-zebra transition 
(left), and using /z c = /x|f of the liquid-zebra transition 
(right). The striking feature is that the distributions re- 
veal a number of peaks. We first discuss the OPD of the 
vapor-zebra transition. Here, the left peak (1) reflects the 
vapor phase, i.e. low colloid density, and high polymer 
density (see corresponding snapshot 1). Although not 
visible in the snapshot, we emphasize that the colloid 
density profile r] c (z) of the pure vapor phase resembles 
that of Fig. [31 i.e. there are (small) density modulations 
along the z-direction. The center peak (2) corresponds to 
a mixed state, where a slab of vapor coexists with a slab 
of zebra phase (snapshot 2). Hence, a vapor-zebra in- 
terface is present, and the corresponding density profile 
rj c (z) will schematically resemble Fig. 0Jb). Note that, 
due to periodic boundaries, the number of vapor-zebra 
interfaces is at least two. The right peak (3) corresponds 
to a pure zebra phase (snapshot 3), with a density profile 
resembling the one shown in Fig. [31 i.e. featuring large 
density oscillations. The meaning of the peaks in the 
OPD of the liquid-zebra transition follows analogously. 
In this case, peak 4 reflects liquid-zebra coexistence, to 
be compared to the profile of Fig. Hfc). Note that the 
density of the zebra phase at the vapor-zebra transition 
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FIG. 14. Analysis of the OPD P(r) c ) obtained at rf v — 1.16, which is above the critical points, but still below the triple point. 
The system size used is L z = 2A = 20, and L x — L y — 12. Shown on the left is P{rj c ) at the vapor-zebra transition, while on 
the right P(rj c ) of the liquid-zebra transition is shown. The salient features are a number of peaks, whose meanings become 
clear upon inspection of snapshots. In the snapshots, the laser field is along the longest edge of the box; colloids are shown as 
black, polymers as white. 



differs from that of the liquid-zebra transition. 

Also of extreme interest are state points "between the 
peaks" in the OPD (Fig. IT5j) . Here, we keep rf p = 1.16, 
but choose larger lateral box extensions, L x — 30 and 
L y = 10, while L z = 2A = 20. In Fig. ITWa). the loga- 
rithm of the OPD is shown, using /i c = /i™ of the vapor- 
zebra transition; note that \nP(r) c ) may be regarded as 
minus the free energy of the system. As in Fig. 1141 three 
peaks are visible: their meaning is the same as before. 
The snapshot of Fig. [T5l b') was taken at i] c =0.1, which 
is between the center and right peak of the OPD. Again, 
vapor-zebra coexistence is observed, but the key differ- 
ence with the coexistence state points of Fig. [14] is that 
one of the periods of the field is only partially filled. 
Hence, in addition to a vapor-zebra interface perpendicu- 
lar to the field, there is a smaller interface parallel to the 
field, indicated by the shaded area A. This is the "hid- 
den" interface, whose presence was already implied by 
the DFT calculation. Note that In P{rj c ) around r\ c = 0.1 
is essentially flat. Hence, once a partially filled slab has 
formed, it can be filled without any cost in free energy. 
This can be understood from the schematic snapshots of 
Fig. ll5f c). which show top-down views (i.e. looking along 
the z-direction) of the partially filled slab; the lateral 
area of the slab equals L x x L y , while the slab thickness 
equals A/2. The snapshots 1,2,3 in (c) correspond to 
state points at the minimum between two peaks in the 
OPD, but with r\ c increasing from left to right (schemati- 
cally resembling the "path" 1 — > 2 — > 3 in Fig.HSTa)). In 
the first snapshot, a droplet of colloidal liquid has con- 
densed. The droplet is cylindrical in shape; note that 
the area of the "hidden" interface in this configuration 
equals the circumference of the circle times the slab thick- 
ness. In the second snapshot, the droplet has grown so 
large it interacts with itself through the periodic bound- 



aries, yielding two slab domains with two interfaces (the 
snapshot of Fig. IToT b) resembles this situation). Since 
L x ^S> L y , the "hidden" interfaces form perpendicular to 
L x ; the shaded region A marks the area of one of them. 
Since the free energy around the minimum of the OPD 
is flat, it follows that the interfaces do not interact^, 
and so we obtain for the surface tension of the "hidden" 
interface^ 

7h = AF/2A, (17) 

with A = ALy/2 (the factor 1/2 in Eq. (JTTJ) is a conse- 
quence of periodic boundaries, which lead to the forma- 
tion of two interfaces). For rf v = 1.16, we obtain 7h ~ 

0. 1fcsT/cr 2 , which significantly exceeds the vapor-zebra 
surface tension. Finally, by increasing rj c even further, 
one obtains the third snapshot, featuring a droplet of col- 
loidal vapor. Note that Fig. IT5l c) is just the "standard" 
droplet condensation transition in a two-dimensional sys- 
tem with periodic boundaries^. 

Having understood the arrangement of the phases in 
the coexistence region, we expect the OPD at the vapor- 
zebra transition to scale with system size conform Fig. 1161 
We assume a L x L x L z periodic box, L z = n\, with L 
large. The two dominating peaks correspond to the pure 
vapor (V) and zebra (Z) phase. The intermediate peaks 

1, 2, . . . , n — 1 correspond to states where vapor and ze- 
bra coexist, with each period of the field completely filled 
with either one of the phases. For each additional period 
of the field, one extra peak arises! The states a, a' , . . . at 
the minima also correspond to vapor-zebra coexistence, 
but where one period of the field is partially filled, im- 
plying the presence of "hidden" interfaces (Fig. [TSTb)). 
In the limit L z — > oo, one thus obtains an infinite se- 
quence of intermediate peaks, separated by "distances" 
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FIG. 15. (a) The OPD obtained for r\ v = 1.16 at the vapor-zebra transition; shown is lnP(i) c ) versus r\ c . The barrier AF 
reflects the interfacial free energy of the "hidden" interface (see details in text), (b) Snapshot taken at rj c — 0.1, which is 
between the center and right peak of the OPD; the laser field propagates along L z . Clearly visible is one completely filled 
slab, and one partially filled slab. The lateral area A of the "hidden" interface is also indicated, where A is the wavelength of 
the field (colloids are shown as black, polymers as white). The schematic snapshots in (c) show the partially filled slab along 
the "path" 1 — !> 2 — !> 3 of (a), where the 2-direction is perpendicular to the plane of the paper (dark regions correspond to 
colloid-rich domains). 



d oc l/L z . The barrier AF reflects the free energy cost 
of the "hidden" interface; AG that of the vapor-zebra 
interface. As L becomes large, we thus expect^ 

AF = 7h A£, AG = 2 7vz L 2 . (18) 

In the thermodynamic limit L — > oc, AG dominates: the 
intermediate peaks then become suppressed, and only the 
pure phase peaks (V,Z) remain. The OPD thus becomes 
bimodal, as it should since the transition is first-order 
between two phases^. The behavior of the OPD at the 
liquid-zebra transition follows analogously, although the 
precise values of 7 h will differ. 

It now becomes clear why 7vz and 7 i z cannot reveal 
critical behavior. The critical behavior was shown to 
be effectively two-dimensional, implying that the singu- 
lar part of the interfacial free energy is due to line ten- 
sion, i.e. proportional to L. As Eq. (|18[) shows, only the 
"hidden" interface reveals this required scaling. Conse- 
quently, 7 h becomes critical, while 7vz , 7 i z do not. 

Note that the OPDs for finite L do not conform to 
Fig. [Tni For instance, in Fig. Q31 the coexistence peaks 
(2,4) exceed those of the pure phases (1,3, 3', 5). This 
is due to the extremely small values of 7vz , 7 i z . For the 
system sizes L accessible in our simulations, AG is essen- 
tially zero, meaning that the intermediate peaks are not 
suppressed. From the finite-size OPD, we thus obtain 
indirect confirmation of the DFT prediction that 7vz , 7 i z 
are extremely small. A second consequence is that the 
tendency of the system to macroscopically phase separate 
is weak. To test this assertion, we performed a semi- 
grand canonical simulation at r\ c = 0.07 and rf p = 1.16, 



using an extremely elongated box with L x — L y = 5, 
L z = 50A. The reader can verify in Fig. [U[b) that this 
state point is deep inside the vapor-zebra coexistence re- 
gion. Consequently, we expect macroscopic phase sep- 
aration, implying the formation of / = 2 vapor-zebra 
interfaces (since the system is periodic). In Fig. [TH1 
we have collected a histogram of observed / values, ob- 
tained in a long simulation run. The key message is that 
the number of interfaces far exceeds two, providing fur- 
ther confirmation that 7vz is small. In some sense, a 
L x L x (L z = nX) system resembles a set of i = 1, . . . , n 
slabs; to each slab we may assign a spin variable, say, 
Si = —1 when the slab is filled with vapor, and Si = +1 
when filled with zebra (further justification of assigning 
spin variables ±1 to slabs follows from the DFT profiles 
of Fig. HJb), which show that the vapor-zebra interface 
is extremely sharp). When two neighboring slabs have 
different spin values, a vapor-zebra interface exists be- 
tween them, which raises the free energy by an amount 
7vz L 2 . This is just the ID Ising models, with Hamil- 
tonian Hi sing ,iD = ~JYa=i s i-i s ii s o = s„, and cou- 
pling constant J — 7vz L 2 /2. In the limit L — > oo, one 
thus recovers the zero-temperature ID Ising model, and 
only here the system will macroscopically phase separate. 
However, due to the small value of 7vz and the finite sys- 
tem size L, it is clear that our simulations are far removed 
from this limit, which also explains the result of Fig. [T5] 
In fact, the solid curve in Fig. IT%1 shows the distribution 
H(I) for the ID Ising model with n = 50 spins and J = 
(with the constraint that the total magnetization X)i=i s i 
is zero). 
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FIG. 16. Sketch of the logarithm of the OPD at the vapor- 
zebra transition, i.e. using ji c — /i^ 2 , in a L x L x L z periodic 
box, L z = nA, and for L large (solid curve). The two dom- 
inating peaks correspond to the pure vapor (V) and zebra 
phase (Z). The intermediate peaks 1, 2, . . . , n — 1 correspond 
to states where vapor and zebra coexist, with each period of 
the field completely filled with either one of the phases. The 
states a, a' , . . . at the minima also correspond to vapor-zebra 
coexistence, but where one period of the field is only partially 
filled, as in Fig JlSf b). The barrier AG reflects the free energy 
cost of the vapor-zebra interface, AF that of the "hidden" in- 
terface; the respective scaling is given by Eq. (jTHJl . 
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FIG. 17. The OPD P(rj c ) at rf p = 1.22 which is close to the 
triple point; system sizes L x = L y = 8, L z = 2 A = 20 are 
used. The peaks corresponding to the pure vapor, zebra, and 
liquid phase are marked (V,Z,L), respectively. 



Finally, we discuss the OPD at the triple point, where 
vapor, liquid, and zebra coexist. In the thermodynamic 
limit, the OPD becomes triple-peaked, each peak corre- 
sponding to one phase. In Fig.[T71 we show the OPD near 
the triple point for a finite system. We indeed observe 
all three phases (V,Z,L) simultaneously, but the coexis- 
tence peaks are still profoundly present. This once more 
confirms the extremely low values of 7 VZ , 7iz , even near 
the triple point, where they are maximal. To describe 
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FIG. 18. Histogram H(I) of the observed number of vapor- 
zebra interfaces I, obtained during a semi-grand canonical 
simulation inside the vapor-zebra coexistence region. The key 
message is that I far exceeds two, showing that the tendency 
for macroscopic phase separation is weak. This is consistent 
with the DFT prediction that 7 VZ is extremely small. The 
dashed curve shows the corresponding histogram for the (ex- 
actly known) ID Ising model at infinite temperature. 



the coexistence in terms of spin variables, as done above, 
now requires 3-state spins, which might induce ID 3-state 
Potts behavior at the triple point (this could be an inter- 
esting topic for further study). Above the triple point, 
only vapor and liquid can coexist, and here the OPD is 
bimodal again2&. In the liquid-vapor coexistence region, 
we could again map the system onto the ID Ising model, 
but with coupling constant J = 7i v L 2 /2. We have veri- 
fied that, due to the substantially larger value of 7i v , the 
tendency of the system to phase separate is now much 
stronger. Of course, for the bulk AO model, the mapping 
onto the ID Ising model does not apply (in this case, 
the external field, Eq. (j4|), which ultimately supplies the 
underlying ID lattice structure, is absent). 

V. CONCLUSIONS 

In conclusion, we have studied fluid phase separation 
inside a static one-dimensional oscillatory external field. 
The actual DFT calculations and simulations were per- 
formed for the Asakura-Oosawa model of colloid-polymer 
mixtures, but we expect that our findings will apply 
to any three-dimensional fluid with a bulk liquid-vapor 
critical point. As was already established in a previous 
work 48 , the external field "splits" the bulk critical point 
into two new critical points, and one triple point. This 
leads to a phase diagram with three coexistence regions, 
featuring (1) vapor-zebra coexistence, (2) liquid-zebra co- 
existence, and (3) liquid-vapor coexistence. All three 
phases (vapor, liquid, zebra) are characterized by den- 
sity modulations along the field direction, but the modu- 
lations are most pronounced in the zebra phase. The im- 
proved DFT calculation of the present work shows that 
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the temperatures of the two critical points differ slightly 
from each other. In addition, we calculated the surface 
tensions associated with all three coexistence regions, and 
found the vapor-zebra and liquid-zebra tensions to be ex- 
tremely small. The DFT calculation also reveals that the 
latter surface tensions do not yield the expected mean- 
field critical exponents (even though our DFT is a mean- 
field theory). 

Computer simulations and finite-size scaling confirm 
all the trends predicted by the DFT. The reason that 
the vapor-zebra and liquid-zebra tensions do not show 
critical behavior is due to the fact that the external 
field divides the system into a series of effectively two- 
dimensional slabs, stacked on top of each other along 
the field direction. The critical correlations diverge only 
in directions perpendicular to the field, and the corre- 
sponding surface tension is one arising from phase coex- 
istence inside single slabs. A surprising finding is that 
the critical behavior is confined to certain slabs only; de- 
pending on the critical point, either the low or high den- 
sity slabs become critical. Along the field direction, and 
above the critical points, the arrangement of slabs can be 
conceived as a one-dimensional Ising chain, at effectively 
zero temperature in the thermodynamic limit, whereby 
each slab represents one Ising spin variable. Hence, there 
will be macroscopic phase separation in all three coex- 
istence regions, but only in the limit where the lateral 
extensions L of the system become large. Regarding the 
vapor-zebra and liquid-zebra coexistence regions, the ten- 
dency to phase separate is particularly weak, due to the 
extremely low vapor-zebra and liquid-zebra surface ten- 
sions. According to our DFT calculations, the latter ten- 
sions are of the order 10 -6 fc^T '/ 'a 2 . This value is too 
low to be quantitatively measured in simulations. How- 
ever, the weak tendency of the system to macroscopically 
phase separate, as observed in our simulations, does con- 
firm that the latter surface tensions must be extremely 
small. 

Our results could be verified in real-space experiments 
of colloid-polymer mixtures using, for instance, confocal 
microscopy^ 7 -. In fact, bulk criticality in these systems 
has already been analyzed in this manner—. The inclu- 
sion of a standing optical field appears to be a feasible 
extensio n 49 ' 69 . In the presence of such a field, the much 
weaker tendency of the system to macroscopically phase 
separate should be easily detectable. 

A remaining puzzle is why our finite-size scaling anal- 
ysis does not reveal two-dimensional Ising critical ex- 
ponents. Of course, the division of the system into 
slabs is not absolute: particles can still diffuse be- 
tween slabs. Perhaps this modifies the universality class, 
but the underlying theoretical mechanism remains yet 
to be elucidated 7 ^. We are currently planning simu- 
lations of the lattice Ising model to address these is- 
sues (the simplicity of the latter model probably allows 
for a more accurate scaling analysis using larger system 
sizes). It would also be interesting to extent the anal- 
ysis to external potentials more complicated than the 



one of Eq. ((4]). Examples include a superposition of sev- 
eral waves resulting in two-dimensional periodi c 51 ' 71 " — 
or quasi-crystalline pattern o 76 ' 77 . The phase diagram of 
a system close to its bulk critical point inside these con- 
fining potentials still needs to be explored. Again, the 
question is whether new critical points arise, and to what 
extent the emerging critical behavior is affected by the 
details of the confining potential. 
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Appendix A: Density functional theory background 

The main variables in our density functional theory 
are the one-body densities Pi(r) of colloids (i = c) and 
polymers (i = p), which describe the microscopic be- 
havior of the system for a given set of parameters (tem- 
perature T and chemical potential pi). Based on the 
existence proofi^ that there is a grand canonical free 
energy functional Q(T, p c , p p , [p c , p p ]) which gets mini- 
mal for the equilibrium density, we use the fundamen- 
tal measure approach^ to approximate this functional. 
The grand canonical free energy functional of a colloid- 
polymer mixture in a three-dimensional system can be 
split as 

fi[p c (r),p p (r)] = J] F-Mr)] + 

i— c,p 

J drp, (r) [T4xt,i 0) - m] + ^c XC [pc (r) , p p (r)] , 

with the external potential V^ x t,i(r) acting on component 
i, keeping in mind a general description where both col- 
loids and polymers are inside an external field. We use 
Eq. ([I]) for the external potential V ex t,c acting on the col- 
loids, and set the external potential acting on the poly- 
cxt.p = 0. In the above, •7 r id[Pi(r)] is the 



mers to zero: V r 



free energy of an ideal gas 

•Fid[Pi(r)] - k B T J drp^r) [in (pi(r)Af) - l] , 

including the (irrelevant) thermal wavelength Aj of the 
particles of species i, an external energy part and the non- 
trivial excess free energy J r oxc [Pc(i"), Pp( r )], which results 
from the interactions of the particles. We approximate 
this excess free energy functional as the integral of a free 
energy density <J>({n£(x)}, {«p(x)}) as 

J- cxc [p c (r),p p (r)] = k B T J dr $(K(x)}, {^(x)}), 
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depending on weighted densities 7iJ,(x) given by the con- 
volution of the actual density profiles with weight func- 
tions 

<( x ) = J drp 4 (rX(x- r). 

The set of weight functions (which are independent of 
the density profiles) is given by 

^« = 0(y-r),*4(r) = *(^-r), 



w\(r) - 

w' 2 (r) 

w*(r) 



^(T- r )'^ = ^(Y 



rr 1 

^2 ~ 3 



with r = r|, the step function 0(r), the Dirac function 
S(r), and the identity matrix 1. The weight functions 
are of different tensorial rank, i.e. scalars w\, w\, w\, Wq, 
vectors w 2 ,w\, and a second rank tensor w 2 . The excess 
free-energy density $ = $i + $2 + $3 is written as 

i— c,p 



= - E 



n 2 n 2 n 2 



— n 2 n 2 ■ n 



c,p 



+ - [njn^n* - Tr (n^n*)] J p«*(n§,ng), 



where Lp l - k (-q c ,jf) = P dv ?™ dr)k Fnr,(i] c , i] p ) denotes the 
derivatives of the OD free energy (3F v(r] c ,r/ p ) = (1 — 
rf — r] p ) ln(l — rj c ) + r] c . We obtain the equilibrium den- 
sity profiles p^(r) by minimizing the Gibbs free energy 
functional, 



Sfl[p c (r),p p (r)} 



0. 



Pi(r)=p< 0) (r) 



This yields the Euler-Lagrange or stationarity equations 
(5$ 



dx 



+ ln(Afo(r)) + /3V e *t,i( r ) - Pfa = 0. 



By inserting the equilibrium profiles into the functional, 
we obtain the grand canonical free energy and can thus 
calculate phase diagrams and interfacial properties. 



